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Abstract 

ov 

The eigenvectors of the Hamiltonian Hn of iV-sites quantum spin chains with elliptic 
exchange are connected with the double Bloch meromorphic solutions of the quantum 
continuous elliptic Calogero-Moser problem. This fact allows one to find the eigenvectors 
via the solutions to the system of highly transcendental equations of Bethe-ansatz type 
^ ■ which is presented in explicit form. 

t> 

It is known [1] that for a one-parameter set Ti^ of linear combinations of iV(iV — l)/2 
elementary transpositions {Pjk}, 7~Cn = jiY^iKj^kPU ~ k)Pjk at arbitrary natural N > 
3, one can construct a variety (3 < m < N) of operators which commute with 

Hn- Being applied to SU(2) spin representations of the permutation group, this proves 
the integrability of ID periodic spin chains with elliptic short-range interaction and the 
Hamiltonian 

<H^= J - £ h(j - k){BjB k — 1), (1) 
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where pjy(x), (n(x) are the Weierstrass functions defined on the torus = C/ZN + Zc<j, 
io = ia, a G R+ is a free parameter. 

The symmetry of two limiting cases of this one-parameter model, i.e. the Bethe lattice 

with nearest-neighbor interaction [2] (a — > 0) and long-range ^ sin j^J exchange [3] 
(a — > 00), is now well understood and regular procedures of finding eigenvectors are 
described in the literature [4-7]. At present, a number of impressive results are known 
for both these models. In particular, they include the additivity of the spectrum under 
proper choice of "rapidity" variables [2,3], the description of underlying symmetry [4,5], 
construction of thermodynamics in the limit N — > 00 [9,10]. However, all that still cannot 
be applied to the general elliptic case. 

In the paper [8], I have shown that there is a remarkable connection between the 
eigenvectors of the Hamiltonian of the above model with M down spins and double Bloch 
meromorphic solutions to the quantum continuous elliptic Calogero-Moser problem at the 
special value of the coupling constant, i.e. the eigenfunctions of the differential operator 

ff = -oEz2+E^-4 (3) 

z 3=1 ux i i+k 

This allows one in principle to find an ansatz for the eigenvectors and even try to describe 
them completely if the solutions to (3) are known. This has been done in the simplest 
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nontrivial case M = 3 in [11], where I have used the result for three-particle elliptic 
Calogero-Moser problem [12]. 

At that time, the explicit form of the eigenfunctions of (3) at M > 3 has not been 
known. The situation has been changed after publishing of seminal paper [13] where 
these eigenfunctions have been obtained in the process of constructing solutions to the 
elliptic Knizhnik-Zamolodchikov-Bernard equations. It has been a main motivation for 
this paper in which I shall describe the complete set of the Bethe-ansatz-type equations 
for the eigenvectors of (1) at arbitrary M < N/2. 

The Hamiltonian (1) commutes with the operator of total spin S = \ ^f=i ®r Then 
the eigenproblem for it is decomposed into the problems in the subspaces formed by the 
common eigenvectors of S 3 and S 2 such that S = S 3 = N/2 - M, < M < [N/2], 



H (s)^(M) >=Em ^(m) > 
The eigenvectors \ip^ M ^ > are written in the usual form 



N 



M 



\^>= ]T ^ M (n 1 ..n M )l[s-\0>, 



(4) 



(5) 



n\..n M 



/3=1 



where |0 >= | ft ••• T> is the ferromagnetic ground state with all spins up and the 
summation is taken over all combinations of integers {n} < N such that Y[^ <u {n^ — n v ) ^ 
0. The substitution of (5) into (4) results in the lattice Schrodinger equation for completely 
symmetric wave function ip M 

N M 

PN{n/3 ~ s)lp M {ni, ..71/3-1, S, 71/3+1, ..1l M ) 

s^nx^.UM /3=1 



+ 



M 



ip M (ni, ~n M ) = 0. 



(6) 



The eigenvalues {Em} are given by 
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where Ci(^) is the Weierstrass zeta function defined on the torus 7\ = C/Z + Zuj. 
The solutions to (6) can be found with the use of the following ansatz for i\) M '- 



<PM( n Pl,~ n PM), 

Pen M 



ip M (ni, ..n M ) 

/ M \ 

(p M \ni, ~n M ) = exp f -i ^p v n v j xSK. -n M ), 



where 



Vv 



Pu - 2nN~ 1 l u , l v e Z, 



(8) 

(9) 
(10) 
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ttm is the group of all permutations {P} of the numbers from 1 to M and Xm i s some 
special solution to the continuum quantum many-particle problem 
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M 



(3=1 
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2 ^ dxl 



E Pn(xp - x x ) 



E M {p) 



Xm(xi,..x m ) = 0. 



(11) 



It is specified up to a normalization factor by the particle pseudomomenta (pi, .-Pm)- The 
standard argumentation of the Floquet-Bloch theory shows that due to perodicity of the 
potential term in (49) Xm obeys the quasiperiodicity conditions 



Xm(xi, -X/3 + N, ..x M ) = exp(ippN)xM (xi, ..x M ), 



(12) 



Xm(^i, -xp + w, ..x M ) = exp(2iriqp(p) + ippUj)x M \xi, -x M ), < Ue{q p ) < 1, (13) 

1<P<M. 

The eigenvalue Em(p) is some symmetric function of (p±, .-Pm)- As will be seen later, the 
set {qp(p)} is also completely determined by {p}. 

It turns out [8] that the equation (6) with the use of (8), (9) can be recast in the form 
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ip {p \n P1 , ..n PM ) = 0, 



where 



//9(p) = 2^(p)Ci(l/2)-Ci(^(p)), 



(14) 

(15) 
(16) 

(17) 



where pi(x), (\{x) are the Weierstrass functions defined on the torus T\ = C/Z + Zou. 

Turning to the definition (9) of (p^ p \ one observes that each term of the left-hand 
side of (14) has the same structure as the left-hand side of the many-particle Schrodinger 
equation (11) and vanishes if Em and fp(p) are chosen as 



fpip) = -iPp, 
E M = E M (p) 



/3 = 1,..M, 

M 



(18) 
(19) 



E^(p)- 

0=1 

One can see from (15-19) that it remains now to find the explicit dependence of {q} 
and Em on {p}. It can be done by using the results given in [13] where the explicit form 
of Xm ( x ) has been indicated. In suitable notations, it reads 



Xm\x) ~ exp(i PP X P) E II °T j ( 



M 



0=1 



; c(s(fc)) x c( 



I(fc))+1 )(*«o-)-*-o-+i)). 



(20) 
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where m = M(M — l)/2, c is non-decreasing function c : {l,..,m} — > {1,..,M — 1} 
such that |c _1 {j}| = M — j, l(s) is an integer which is defined for the permutation s by 
the relation x cis{1 )) +1 d/dx c ( sW )...x c{s{m )) +1 d/dx cis{m ))X^ = l(s)(xi...x M ), {£} is a set of m 
complex parameters obeying m relations [13] 

2 P(tj -t{)-2 ]T P( f j - f i) + M5 cj,iP{ t j) = KPcU) - Pc(j)+i), (21) 

l:\c(l)-c(j)\=l hl^j, c{l)=c{j) 

Pit) = CatW - ^(N(N/2)t, 

and 

a w (t)=exp((2/N){ N (N/2)wt) £7 ' v( "' "° 



0-Jv(w)0-7v(t)' 

<tat being the Weierstrass sigma function on TV The elementary building blocks of the 
X function obey the useful quasiperiodicity relations 

flWiv(t) = cr«,(<), ^+o;(i) = e 2 ^ /Jv a w (t). (22) 

One can see that in this construction the color function c(j) is of crucial role. It is useful 
to write it explicitly. Namely, define for every k=l,..M — 1 the segment Sk 

(fc-l)(2M-fc) + x < J < fc(2M-fc-l) _ (23) 

Then some calculation shows that 

c(j)=k if J6S, (24) 

The main advantage of the explicit form of x function is that it allows to find the 
second set of relations between the Bloch factors {p}, {q}. It is easy to see from (21) that 
{p}'s in the definitions (12) and (20) are the same. The problem consists in calculation 
of {q}. To do this, it is not necessary to analyze each term in the sum over permutations 
in (20) since all of them must have the same Bloch factors. It is convinient to choose the 
term which corresponds to the permutation 

so: s (j)=m+l-j, j = l,..m. 

After some algebra, one finds that this permutation gives nontrivial contribution to the 
sum (20) with 1(sq) = M\(M — 1)!...2!. Moreover, with the use of explicit form of the 
color function (23-24) one finds 

c(s (l)) = M-q if q(q- l)/2 + 1 < I < q(q + l)/2. 

Now the problem of calculation of the second Bloch factors reduces, due to second re- 
lation (21), to some long and tedious, but in fact simple calculations of the product of 
factors which various a functions acquire under changing arguments of x function to the 
quasiperiod u. The final result is surprisingly simple, 

qp{p)=N~ 1 [ £ U- Yl *i)> K/3<M-1, (25) 

\l:c(l)=P l:c(l)=P-l J 



4 



with the first and second term being omitted for (3 = M and (3 — 1. 

The equations (25), together with (18) and (21), form a closed set for finding Bloch 
factors {p}, {q} at given integers {Ip} G Z/MZ and determining the eigenvalues of the spin 
Hamiltonian (1,2) completely. The corresponding eigenvalue of the continuum M-particle 
operator (11) is given by [13] 



Em(p) = 



2M(M 



N 



M 

/j=i 



5Z(2<J c (fc), c (j)-P'(*fc 

k<l 



U) - &\c{k)-c(l)\,lF{t k - ti)) 



m y: 

c(k) = 



F(t k ) 



(26) 



where 



F(t) = -p N (t) + (( N (t) - 2/N( N (N/2)) 2 + 4/N( N (N/2). 



This allows one to find, via (7) and (19), the explicit form of the eigenvalues of spin 
Hamiltonian (1,2). It is worth noting that for their real calculation one has to solve the 
Bethe-type equations (18), (21), (25) at first. 

In conclusion, it is demonstrated that the procedure of the exact diagonalization of 
the lattice Hamiltonian with the non-nearest-neighbor elliptic exchange can be reduced 
in each sector of the Hilbert space with given magnetization to the construction of the 
special double quasiperiodic eigenfunctions of the many-particle Calogero-Moser problem 
on a continuous line. The equations of the Bethe-ansatz form appear very naturally as 
a set of restrictions to the particle pseudomomenta. The proof of this correspodence 
between lattice and continuum integrable models is based only on analytic properties of 
the eigenfunctions. One can expect that the set of spin lattice states constructed by this 
way is complete. This is supported by exact analytic proof in the two-magnon case. 

The analysis of explicit form of the equations (21) available for M = 2, 3 shows that 
the spectrum of the lattice Hamiltonian with the exchange (1) is not additive being given 
in terms of pseudomomenta {p} or phases which parametrize the sets {p, q} [11]. For 
arbitrary M, this can be seen directly from (26). The problem of finding appropriate set 
of parameters which gives the "separation" of the spectrum remains open. It would be 
also of interest to consider various limits (N — > oo, a — > 0, oo) so as to recover the results 
of the papers [2,3] and prove the validity of the approximate methods of asymptotic Bethe 
ansatz. 



5 



References 

1. V.I. Inozemtsev. Lett. Math. Phys. 36,55 (1996) 

2. H. Bethe. Z.Phys. 71,205 (1932) 

3. F.D.M. Haldane. Phys. Rev. Lett. 60,635 (1988); B.S. Shastry, ibid 60,639 (1988) 

4. L.D. Faddeev, in Integrable models of 1+1 dimensional quantum field theory, Else- 
vier, Amsterdam, 1984 

5. D. Bernard, M. Gaudin, F.D.M. Haldane and V. Pasquier. J.Phys. A26,5219 (1993) 

6. M. Fowler and J. Minahan. Phys. Rev. Lett. 70,2325 (1993) 

7. B.S. Shastry and B. Sutherland. Phys. Rev. Lett. 71,5 (1993) 

8. VI. Inozemtsev. J.Phys. A28,L439 (1995) 

9. F.D.M. Haldane. Phys. Rev. Lett. 66,1529 (1991) 

10. M. Takahashi. Progr. Theor. Phys. 46,401 (1971) 

11. V.I. Inozemtsev. J.Math. Phys. 37,147 (1996) 

12. J. Dittrich and V.I. Inozemtsev. J. Phys.A: Math. Gen. 26,L753 (1993) 

13. G. Felder and A. Varchenko. Int. Mat. Res. Notices 5,222 (1995) 



6 



